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ABSTRACT 

We use Doppler imaging techniques to determine the dependence of starspot rotation 
rates on latitude in an homogeneous sample of young, rapidly-rotating solar analogues. A 
solar-like differential rotation law is used, where the rotation depends on sin 2 (6), where 9 is 
the stellar latitude. By including this term in the image reconstruction process, using starspots 
as tracers, we are able to determine the magnitude of the shear over more than one rotation 
cycle. We also consider results from matched filter starspot tracking techniques, where indi- 
vidual starspot rotation rates are determined. In addition we have re-analysed published results 
and present a new measurement for the K3 dwarf, Speedy Mic. 

A total of 10 stars of spectral type G2 - M2 are considered. We find a trend towards 
decreasing surface differential rotation with decreasing effective temperature. The implied 
approach to solid body rotation with increasing relative convection zone depth implies that 
the dynamo mechanism operating in low-mass stars may be substantially different from that 
in the Sun. 

Key words: Line: profiles - Methods: data analysis - Techniques: miscellaneous - Stars: 
activity - Stars: atmospheres - Stars: late-type 



1 INTRODUCTION 

ICarringtonl 1 186(1) was the first to attribute the scatter in rotation 
rates of sunspots to a differential rotation of the solar sur face. It 
would however be a further 50 years before iHald ^908) would 
discover the existence of magnetic fields in Sunspots. Helioseis- 
mological studies of acoustic-wave propagation through the solar 
interior confirm that the equatorial regions rotate more rapid ly than 
the po lar regions throughout the outer convective zone (e.g. lSchod 
1998). Radial shear is most pronounced at the tachocline which is 
located at the interface between the convection zone and radiative 
solar interior. 

The tachocline provides a relatively stable environment for 
the amplification of magnetic field, being located within the re- 
gion where convective cells can sink and overshoot into the solar 
core. In this region of strong radial shear, the fi-effect stretches 
poloidal (North-South field) field into toroidal (East-West) field. 
Toroidal field components entrained in convecting fluid elements 
rotate under the influence of Coriolis forces as they rise and fall, 
regenerating th e poloid al field via the a- effect <Parkedll955l 
iKrause & Raedledll98Cl) . One goal of solar dynamo theory is to 



combine these effects into a self-consistent physical model that can 
reproduce the 1 1-year solar activity cycle in detail. 

Magnetic activity in stellar counterparts of the Sun gives rise 
to periodic optical variability as starspots rotate into and out of 
view. This variability is also seen in optical and UV emission 
lines from the chromosphere, and in the coronal X-ray luminos- 
ity. All these indicators are correlated with photospheric magnetic 
field strength on the Sun, and increase with rotation rate in stars of 
a given mass. 

Long term photometric monitori ng of stars pot modul ation in 
magne tically-active binary stars by lHalll ll99ll) and iHenrv et alJ 
1 1995) (H95) has revealed secular changes in rotation period which 
are strongly indicative of rotational shearing. As a stellar activity 
cycle progresses, starspots should emerge at different latitudes. If 
the star rotates differentially, these changes in latitude will produce 
long-term scatter in the photometric rotation periods measured at 
single epochs. The correlation between this scatter and the mean 
rotation period indicates that AS1 oc ^ - 24 ± - 06 (where O, is angu- 
lar velocity and AQ the corresponding scatter in radians day -1 ), 
and thus that AQ is almost independent of rotation (H95). 

The Mount Wilson Survey of chromospheric Ca II H & K 
emission-line fluxes has shown that many solar-like stars exhibit 



2 J.R. Barnes, A. Collier Cameron, J.-F. Donati, D.J. James, S.C Marsden, P. Petit 



Differential rotation x 2 landscape for Speedy Mic (2002 Jul/Aug) 

0.05 | 1 1 1 1 , , 1 . 1 1 ^ , 1 — n 1 1 1 1 , 1 1 1 1 1 



0,04 - 




0.02 - 



0.01 ' ' 1 1 1 ' 1 1 1 1 ^ — ' ' ' — J ' ' ' 1 ' 1 1 1 1 — 

16.534 16.535 16,536 16.537 

E! [rad.day~ 1 ] 

Figure 1. The figure here shows a x 2 landscape for the 2002 July 18, 19 
& 23 Speedy Mic (HD 197890) timeseries, revealing that rotation rate and 
differential rotation are correlated. The minimum point (marked by a cross) 
gives the best fitting differential rotation and related errors through the 1- 
parameter l-<r confidence interval (inner contour). Also shown are the 2- 
parameter 1-rr and 2.6-rr (99%) confidence intervals. 



multi-periodic chromospheric emission flux variations (Wilson 
ll978tlBaliunas et alll 995). The short-period modulation (of order 
days to weeks) gives the rotation period of the dominant active lat- 
itude at any epoch. The longer-period modulations correspond to 
the solar activity cycle. Over timescales of many years, the rotation 
period is itself found to show a characteristic sinusoidal variation 
in many stars of spectral type F to K. Cyclic changes in the period 
are attributed to a st ellar dynamo driven by differential rotation. 
iDonahue et all < 19961) (DBS 96) have shown that AS1 oc ft 7±0 \ a 
result which is not consistent with the findings of H95, even within 
the quoted uncertainties. The sample of H95 comprises mostly bi- 
nary systems, although a correction term for filling factor is in- 
cluded in their analysis. It should be emphasised here that no lat- 
itudinal information about the starspot groups or active regions is 
assumed in the Afi /Q. proportionality. A more significant explana- 
tion for the discrepancy of the results presented by H95 and DBS96 
is related to this issue and will be discussed in fJ5| 



2 METHOD & RESULTS 

The work reported here employs the Doppler imaging technique, 
which offers a more direct approach than the above methods to 
measuring stellar surface rotation rate as a function of latitude. This 
sheared image method allows us to make an indirect image of the 
stellar surface using a model which include s the effects of d iffer- 
ential rotation and has been investigated bv lPetitetalJl2002t) . The 
image-data transformation incorporates a simplified solar-like dif- 
ferential rotation law, where rotation is proportional to the square 
of the sine of the stellar latitude, 0, i.e. = Q eq — Afi.sin 2 (f?). 
We determine the value Af2 of the stellar surface shear that gives 
the best fit to the data secured over several stellar rotations. Since 
spots are present at a variety of latitudes, their individual rotation 
rates trace the differential rotation pattern. Fig.Qdemonstrates how 
the goodness of the final fit to the data is used to determine the equa- 
torial rotation rate. Here we present the results from data taken with 
the Anglo Australian Telescope and the University College London 
Echelle Spectrograph. Observations of the 0.38 d period K3 dwarf 
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calculated Af2 using the sheared image method (see text for details). 
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Speedy Mic (HD 197890) were secured on the nights of 2002 July 
18, 19 & 23 and will be presented in detail in a separate publication. 

Additionally, Collier Cameron & Donati (2002) applied a 
matched filter analysis technique, allowing the rotation rates of in- 
dividual starspots to be determined directly from the timeseries 
spectra. We include those measurements for AB Dor in addition to 
results from the sheared image method (D onati et all2 003). 

This technique has been applied to an homogeneous sam- 
ple of ten young, rapidly rotating stars with spectral types rang- 
ing from G2V to M2V. The measured shear, B-V colours and 
referenc es are listed i n Table Q The model-atmosphere relation 
given bv lBessellet~aT]ll998l) was used to obtain the effective sur- 
face temperature from the observed B-V colours. In two cases, 
we have re-determin ed AO us ing the sheared image method. 
For PZ Tel iBarnes et alj |2000) we initially used image cross- 
correlation which yielded a slightly higher value for the lati- 
tude depen dent shear. We have a lso re-analysed the results from 
IBarnes & Collier Cameror] fcOOll) for the M1-2V EY Dra (RE 
1816+541). In this latter instance, we initially used splines to in- 
terpolate the grid search of AQ vs Q, which gave a very non- 
uniform x 2 landscape with more than one minimum. We now fit 
a quadratic function in both dimensions eliminating local fluctua- 
tions and yielding a single global minimum. The combined results 
from all studies are plotted in Fig. [3] and reveal a marked trend 
towards lower differential rotation with decreasing stellar effective 
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Figure 2. Differential rotation (Af2) as a function of stellar temperature. The circled dot represents the Sun. Independent measurements, made in different 
years, are shown for three of the stars in our sample: AB Dor, LQ Hya and R58. A dashed line connects pairs of points for LQ Hya and R58. 



temperature. We find a power-law dependence of the surface shear 
on stellar surface temperature: 



AO oc T 



(1) 



Those stars with more than one measurement of the surface 
shear show that the magnitude of the scatter is greater than that of 
the unc ertainty on a si ngle measurement, indicating secular vari- 
ability ICollier Cameron & Donatl2002tlDonati et all2003l) in ei- 
ther the differential pattern itself or in the depths at which the fields 
giving rise to the tracers are anchored iDonati et al. 2003). Never- 
theless LQ Hya poses the most problematic point in the measured 
relation. Careful analysis suggests that the measured shear values 
are reliable. We note that LQ Hya is the slowest rotator in the sam- 
ple, yielding relatively poor spatial resolution (~16°) at the stellar 
equator. Short-term evolution of unresolved spot groups may lead 
to a systematic mis-determination of the differential rotation. We 
discuss this issue and secular variations of AQ further in [J3] 

Fig. [3] shows the dependence of differential rotation, AQ, 
on equatorial rotation rate, Q, for the stars in our sample, along 
wi th a vailable d ata f r om ot her studies. We discuss the results 
of iReiners & Schmittl <2003l) (RS03) in §3\ The unweighted least 
squares fit (i.e. neglecting errors or limits) yields a slope indicating 



that AD. oc Q" 



The values for AQ derived from Doppler 



imaging studies and the solar value (circled dot) represent the full 
difference in rotational frequency between the pole and the equa- 
tor, while the values for AQ derived from DSB96 and the H95 are 
lower limits. In other words, they represent the spread in rotation 
rates only over the range of latitudes where active regions erupt dur- 
ing the stellar cycle. The scatter in AQ at any given rotation rate is 
significantly greater than the trend in mean AQ over the range of 
Q plotted. 

Despite the limitations of an unweighted fit, the power law 



relationship between AQ and Q is consistent with the findings of 
H95. Here we can see that the slope in the DSB96 data alone is 
greater. A possible explanation is that there is a spectral type de- 
pendency on the range of latitudes at which magnetic flux errupts. 
The faster rotators in the DSB96 sample tend to be later spectral 
type, while the slower rotators are of earlier spectral type. The lat- 
ter results are in agreement with RS03 who obtained significant 
measures of differential rotation for early G dwarfs and F dwarfs 
only. Despite the high degree of scatter, there is some evidence for 
spectral type dependence in Fig.|3|when we compare the results of 
RS03 (high AQ) with those presented for the M dwarf stars using 
Doppler imaging studies (low AQ). It is this dependence which we 
have discerned for the first time in Fig.|2|for a narrow range of Q. 



3 DISCUSSION 

RS03 found recently that main-sequence stars of spectral types F 
and G exhibit significantly greater differential rotation than we find 
in the sample of lower-mass stars presented here. The high rota- 
tional shear of the F and G stars manifests itself in the Fourier trans- 
forms of their rotationally broadened line profiles, allowing precise 
measurements of shear from line profiles alone for stars with inter- 
mediate rotation rates. We have discussed the enhanced rotational 
shear (Fig.|5J in the thi n con v ective zones of these stars in ij2] 

Kitchatinov & Riidiaet 1 1999) have modelled internal fluid 
motions to obtain the surface differential rotation in various stel- 
lar types. Their model predicts that the surface differential rota- 
tion rate should vary only weakly with axial rotation period and 
spectral type in the G2 - K5 range. A K5 star is predicted to show 
approximately 2/3 the surface shear of a G2 star. This is in quali- 
tative if not quantitative agreement with our finding that differen- 
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Figure 3. Differential rotation (logl()(Af2)) as a function of loglO(C). Results from various studies: Doppler imaging (see Table 1), DSB96, H95 and RS03. 
The Sun is represented by the circled dot. The dashed line is an unweighted fit to all points with slope 0.15 ± 0.10. 



tial rotation is an order of magnitude lower at K5 than at G2. The 
stars in our sample rotate more rapidly that the fastest stars in the 
models of Kitchatinov & Riidiger, and it is not clear what effects 
such rapid rotation would impose on model solutions. These au- 
thors however found that AO, <x O on average for a G2 star 
and AO, cx f2 -0 05 for a K5 star, in reasonable agreement with the 
unweighted fit (Afi cx n 015±0 10 ) shown in Fig.|3l 

That scatter is seen for some stars with more than a single 
measurement requires further comment. Being the brightest mem- 
ber of its class, AB Dor is the most studied single rapid rota- 
tor. While there is a slight offset between AS7 as measured us- 
ing starspot tracking and the sheared image techniques, the over- 
all seasonal variations show that the equator-pole lap time varies 
from 70 d to 140 d. An analogue of the torsio nal oscillations seen 
at the stellar surface l Howard & Labonte 1980) and within the con- 
vection zone (Kosovichev & Schou 1997) may be responsible, but 
the variation in shear is between one and two orders of magnitude 
greater than the solar variations which amount to 5 m s" 1 dev iations 
from fixed differential rotation. It has been shown t Apoleaate Il992t) 
that the transport of angular momentum may change in response to 
changes in the magnetic flux threading the convection zone. This 
behaviour may be expected to be periodic if magne tic activity cy- 
cles are present. ICollier Cameron & Donatil 120021) discuss these 
possibilities and show that the observations of AB Dor are consis- 
tent with the mec hanism put forward by Applegate. Alternatively, 
iDonati et alj J2003I) have suggested that the starspots which trace 
differential rotation may be anchored at different depths in the 
convection zone. This implies that the magnetic field is generated 
through the action of a distributed dynamo in rapid rotators. 

In the case of the more slowly rotating LQ Hya however, 
we see more than an order of magnitude variation in Afi. We 
have briefly suggested (f|2j a possible cause of mis-determination 



of differential rotation shear in more slowly rotating stars, but 
note this argument does not apply to the rapidly rotating R58 
(v sini= 92 km s -1 ) which shows considerable variation in shear. 
Clearly multi-epoch measurements of differential rotation are re- 
quired for each star if this phenomenon is to be investigated fur- 
ther. Our confidence in the Afi vs temperature trend is affirmed 
by the observation that a statistical sample of both M dwarfs and 
mid-late K dwarfs all ex hibit a small AS! while the findings of 
Reiners & Schmitt 1 2003) for F dwarfs indicate greater Afl than 
for those stars in our sample. 

The decrease in differential rotation that we observe toward 
later spectral types suggests that the importance of the a-effect rel- 
ative to the Sl-effect should increase strongly toward lower-mass 
stars. The a-effect should also increase with increasing rotation 
rate, since it arises from the action o f Coriolis forces on a con- 
verting fluid iKrause & RaedleJT9 80?). The rotational shear, how- 
ever, appears to be weakly dependent on rotation rate, and indeed 
decreases strongly towards lower stellar masses. We conclude that 
the a-effect is therefore likely to dominate stellar magnetic-field 
generation in low-mass stars, and in rapidly rotating young stars 
with an even wider range of stellar masses (e.g. Kiiker & Riidiger 
1999) 

Whether a change in the dynamo mechanism occurs as the 
a-effect becomes dominant remains unclear. With no shear bound- 
ary, the nature of the dynamo must change in fully convective stars 
so that an a 2 mechanism with a dynamo distributed throughout the 
convection zone would then be expected to generate and maintain 
magnetic activity (Kiiker & Riidiaer 1999). Nevertheless, the mass 
at which stars become fully convective is unclear, since a change 
in the nature of dynamo is reasonably expected to lead to a break 
in activity-rotation indi cators. Only recently ha s such a disconti- 
nuity been discovered l M ohantv & Ba sri 2003), albeit at a much 
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later spectral type (M9) than stellar structure models normally pre- 
dict. These authors suggest that an increasingly neutral atmosphere 
could res ult in a d rop of magnetic activity at this spectral type. Ad- 
ditionally, Mullan & MacDonald (2001) have shown that the inter- 
nal strength of the magnetic field can modify stellar structure sig- 
nificantly such that stars become fully convective at later spectral 
type than is predicted by standard evolution models. 

While there is no evidence for a discontinuity in activity in- 
dicators at early M spectral types, other observations may be in- 
dicative of a change in the nature of dynamo activity. Although 
coronal X-rays tend to be saturated in the rapidly rotating regime, 
decreasing optical lightcurve amplitude for late K through to early- 
mid-M (Messina et al. 2003) could be suggestive of uniformly dis- 
tributed spots. This may be explained by distributed dynamo ac- 
tivity. Also, the model of ( Kiiker & Rtidieeil ll999l) indicates that 
the magnetic field is steady, and has a bipolar geometry, but with 
the axis inclined at 90° to the rotation axis. While magnetic Zee- 
man imaging may be required to test the prediction of an in- 
clined field, uniformly distributed spots and a steady state are in- 
dicat ed by si milar distributi ons of star spots on the dM1.5 star HK 
Aqr jBarnes & Collier Cameror]|200ll : lBarnes et all2004 at three 
epochs over a period of 1 1 years. 
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